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Abstract. The magnetization process of an S = | antiferromagnet on the kagome 
lattice, [Cu3(titmb)2(OCOCH 3 )6]-H 2 0{titmb= l,3,5-tris(imidazol-l-ylmethyl)-2,4,6 trimethyl- 
benzene} has been measured at very low temperatures in both pulsed and steady fields. We 
have found a new dynamical behavior in the magnetization process. A plateau at one third 
of the saturation magnetization M B appears in the pulsed field experiments for intermediate 
sweep rates of the magnetic field and disappears in the steady field experiments. A theoretical 
analysis using exact diagonalization yields, Ji = — 19 ± 2K and J2 = 6 ± IK, for the nearest 
neighbor and second nearest neighbor interactions, respectively. This set of exchange parame- 
ters explains the very low saturation field and the absence of the plateau in the thermodynamic 
equilibrium as well as the two-peak feature in the magnetic heat capacity observed by Honda 
et al. [ Z. Honda et al, J. Phys.: Condens. Matter 14, L625 (2002)]. Supported by numerical 
results we argue that a dynamical order by disorder phenomenon could explain the transient 
appearance of the M s /3 plateau in pulsed field experiments. 



The study of quantum phenomena in magnetic systems has been a central subject in 
condensed matter research in recent years. Lower dimensional antiferromagnets show many 
interesting phenomena originating from quantum fluctuations. One of the examples of these 
phenomena is the appearance of a singlet ground state and a spin gap to the excited triplet in a 
spin, ,5=1 one-dimensional Heisenberg antiferromagnet (ID HAF) [1], a spin-Peierls system [2] 
and an S=| two- leg spin ladder [3]. Another interesting example is the appearance of a mag- 
netization plateau, a phenomenon in which the magnetization, M, of an antiferromagnet stays 
constant in a finite range of applied magnetic fields, B, before reaching the saturated phase. 
In the last few years, magnetization plateaus have been observed in several quantum antifer- 
romagnets including an S=l ID antiferromagnet with bond alternation [4], the coupled dimer 
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compound NH4CUCI3 [5] and SrCu 2 (B03) 2 [6-11]. The origin of the magnetization plateau in 
these quantum antiferromagnets differs for different compounds since a magnetization plateau 
can be expected to occur from various combinations of dimensionality and frustration of in- 
teractions [12]. Theory predicted that a quantum antiferromagnet on the kagome lattice 
with a nearest neighbor Heisenberg interaction, which is a typical example of frustrated mag- 
nets, would also show a magnetization plateau at one third of the saturation magnetization 
M s [12-14]. There has been, however, no experimental finding of the magnetization plateau 
in kagome antiferromagnets. In this paper, we report the first observation of a magnetization 
plateau in the S=^ kagome antiferromagnet, [Cu3(titmb) 2 (OCOCH 3 ) 6 ]-H 2 0, {titmb=l,3,5- 
tris(imidazol-l-ylmethyl)-2,4,6 trimethylbenzene} (hereafter, abbreviated to Cu-titmb). 

First, we summarize the crystal and magnetic properties of Cu-titmb. The crystal structure 
of this material has been determined by Liu et al. [15] and is shown schematically in Fig. 1. 
As is seen from Fig. 1(a), Cu 2+ ions are placed at the vertices of the triangles. The distance 
between the Cu 2+ ions on a triangle is 7.77 A. These triangles are attached to a hexagon. A 
two dimensional array of such pattern is called the kagome lattice. The kagome layers are 
stacked along the c axis as shown in Fig. 1(b) with the inter layer distance of 7.05 A. We 
define the nearest neighbor exchange interaction, J\, between Cu 2+ spins in a kagome layer 
as the one along the edges of a triangle through H 2 and/or CH 3 COO~. The former path 
is expected to dominate because the latter path is connected by a weak hydrogen bond. The 
second neighbor interaction, J 2 , is defined as the one between a spin and the spins on the 
second neighboring vertices of the hexagon. Since no apparent exchange paths exist between 
spins on the adjacent kagome layers along the c axis, the interaction between layers is expected 
to be much smaller than Ji and J 2 . 

Heat capacity measurements were performed on a powder sample of Cu-titmb [16]. The 
magnetic heat capacity, C m , obtained after subtraction of the lattice heat capacity is shown 
in the inset of Fig. 2 as a function of temperature, T. We see a two peak structure which 
is a feature predicted for an S — \ aniferromagnet on the kagome lattice with the nearest 
neighbor Heisenberg interaction [17-19] and for other frustrated quantum magnets [12]. 

Powder samples of Cu-titmb were grown by the method described in [15]. The material, 
l,3,5-tris(imidazol-l-ylmethyl)-2,4,6-trimethylbenzene, was purchased from the Wako Pure 
Chemical Industries, Ltd. The magnetization measurements were performed at Osaka Uni- 
versity under pulsed high fields in conjunction with a dilution refrigerator. Pulsed magnetic 
fields are generated by discharging the energy stored in a capacitor bank through the coil. 
The duration time of the pulsed field is given roughly by (LC) 2 , where L is the inductance of 
the coil and C is the capacitance of the bank. Since the magnetic field changes sinusoidally 
with time, the field sweep rate, dB/dt, differs for different field values. In the following, we 
denote the field sweep rate as the one obtained at around 1.5 T. To change the sweep rate, 
we have to change either, C, L and/or the peak field. In this experiment, we used two sets of 
capacitor banks and were able to measure the magnetization over a wide range of sweep rates. 
However, we were not able to access the sweep rates below 0.7 T/msec and those between 1.9 
T/msec and 14 T/msec from the technical reason described above. The powder sample was 
packed into a Teflon tube and immersed in the mixing chamber of the dilution refrigerator. 
In the pulsed field measurements, we observe the derivative of magnetization with respect to 
time, dM/dt. Simultaneously, we measure the time evolution of B. Combining these two sets 
of data, we get M vs B curve. Figure 2 shows the magnetization as a function of applied mag- 
netic field measured at 100 mK in increasing and decreasing fields. We see a magnetization 
plateau at M s /3 for 0.5 T < B < 2.5 T in the increasing field run. Then the magnetization 
increases almost linearly with increasing field up to saturation. The magnetization becomes 
saturated above about 10 T with the saturation moment of about 1 fj,B per Cu atom. This 
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value of M s is consistent with the value expected for a Cu 2+ with S — \ and g = 2.20 [15]. 
Surprisingly, no plateau has been observed in decreasing field. Figure 3 shows the magnetiza- 
tion curve measured in increasing fields with different sweep rates. As the field is swept slower, 
the M s /3 plateau becomes less pronounced. When B is swept much faster, we do not see the 
plateau either. The results presented in Figs. 2 and 3 suggest that the M s /3 plateau might 
be metastable. In order to get a magnetization curve in equilibrium, we have measured the 
magnetization of Cu-titmb by the Faraday method in steady fields using a superconducting 
magnet [20] . A powder sample of Cu-titmb was mixed with Apiezon-N grease and put on one 
of the plates of the capacitor. In order to ensure a good thermometry, Cu foils were used to 
link the sample and the mixing chamber. Figure 4 shows the results. In the thermodynamic 
limit, the magnetization increases abruptly with increasing fields from almost zero and no 
plateau is seen. Saturation of the magnetization is observed at B sat ~ 0.2 T. 

In the following, we discuss theoretically the experimental results reported above. The 
proximity of a ferromagnetic instability due to a competition between ferromagnetic and 
antifcrromagnetic exchange interactions is the simplest explanation for the anomalously low 
saturation field, S sat , measured in the static experiments. Such competing interactions may 
be described by a J\ — J2 model on the kagome lattice, with J\ x J 2 < 0. The model reads: 

<i,j> «i,k» 

where the first and second sums run on the first and second neighbors respectively. For pure 
antifcrromagnetic J\ (or J2) couplings and in a finite range of parameters around these sit- 
uations the system is in the special spin liquid phase which has been the subject of many 
studies [17, 18,21]. A large enough second neighbor ferromagnetic interaction drives the sys- 
tem towards the so-called V3 x \/3 three-sublattice Neel order, whereas a combination of 
antifcrromagnetic J\ and J2 induces a q — three-sublattice Neel order [22] . A ferromagnetic 
phase exists over a large range of J\ and J 2 , which become unstable when the antifcrro- 
magnetic second-neighbor interaction becomes large (J2/J1 > —0.35). As will be explained 
below, we found that Cu-titmb is in this range of parameters very close to the ferromagnetic 
phase. For these coupling parameters a numerical calculation up to 36 sites gives the quantum 
eigen-spectra which show a 12 sublattice Neel order with a non coplanar order parameter, in 
agreement with the classical solution. A spin wave calculation remains to be done to give a 
plausible estimate of the reduction of the order parameter by long-wave-length quantum fluc- 
tuations [23]. In the experiments _B sat is very low compared to T max , the temperature at which 
C m shows a broad peak (see Fig. 2): the dimensionless parameter, 9 = (<7/UB-B S at) / (felmax) 
is approximately 0.02, whereas this ratio is expected to be of the order of one in an un- 
frustrated antiferromagnet and is 4.5 in the pure antifcrromagnetic model on the kagome 
lattice [12-14, 17]. In the J\ — J 2 model, 6 decreases when approaching the boundary of the 
ferromagnetic phase [23] . A value of 9 consistent with the experimental one is obtained in the 
vicinity of the ferromagnetic instability either with ferromagnetic J\ and antifcrromagnetic 
J2 with J2/J1 ~ —0.3, or with antiferromagnetic J\ and ferromagnetic J2 with J2/J1 ~ —50. 
Mizuno et al. [24] have calculated the exchange interactions in copper oxides. They showed 
that the exchange interaction between Cu 2+ spins through an oxygen atom is ferromagnetic 
when the bond angle is around 90°. We, therefore, adopt a ferromagnetic J\. 

A reasonable agreement between the experimental data for T max and B sat is obtained with 
J\ = — 19 ± 2K and J2 = 6 ± 2K. This set of J\ and J 2 also predicts a double peak feature 
in C m as observed experimentally (Fig. 2). At the thermodynamic limit, this model has no 
magnetization plateau as will be shown below. 
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Such a competition of interactions induces a very high density of low lying magnetic states. 
Numerous magnetic excitations are thus at very low energy and very small perturbations can 
change the magnetic behavior. We have observed in our numerical study a general phe- 
nomenon that might be relevant to the real physical behavior of the compound. Samples 
with N~ 21 and 27 which do not accommodate the wave vectors at (0, ir), (it, 0) or (it, it) 
in their Brillouin zone frustrate the thermodynamic ground state and display a I plateau as 
shown in Fig. 5. In this figure, is shown the dependence of the field range, B 2 — B\, at which 
the magnetization plateau at m(= M/M s ) = ^ appears on the energy per spin, e(m), in 
this state. At the thermodynamic limit shown schematically in the inset, a cusp in the e(m) 
curve at m=m is the signature of a m magnetization plateau bounded by the two fields Bi 
and B 2 which measure the slopes of the e(m) curve on the left and right of mo [12]. Exact 
diagonalization on finite size samples gives the energy per spin versus magnetization for the 
discrete values m=2S to t/N. These energies are then fitted to two polynomial forms allowing an 
eventual discontinuity of the slopes at the specific point mo . The best fits are obtained with 
a discontinuity at m = | and give the width B 2 — B\ of the magnetization plateau reported 
in the figure. The error bars represent the uncertainties in the parameters of each fit. The 
results for the sizes multiples of 6 (N=18, 24, 30 and 36) exhibit zero width (i.e. no plateau) 
in the thermodynamic limit. The N=21 and 27 samples, which do not allow the wave vectors 
(0,7r), (71-, 0), (7r,7r) have a higher energy per spin and exhibit a plateau at m = |. At any 
magnetization the energy per spin of these samples at T=Q is indeed slightly higher than the 
true thermodynamic ground-state and thus the system is in a metastable configuration. At 
M s /3, this metastable state is associated with an up-up-down (uud) configuration of spins, 
where spins along the horizontal lines in Fig. 1(a) are ferromagnetically aligned in the up 
direction, whereas the extra spins of the triangles are down [12]. Three such configurations 
are possible and thus we have three domains. 

Among all the configurations of spins that are at very low energy at M s /3, purely geometric 
arguments imply that the three uud configurations have in phase space the largest attraction 
basin for fluctuations, either thermal or quantum, and thus a lock-in of the system in such 
configurations can be driven by a large number of very small perturbations. This is a transient 
analogue of the phenomenon of order by disorder [25-27]. In such a scenario there exist in 
the free energy landscape three large metastable wells associated with the three possible uud 
configurations. Two characteristic times of the out of equilibrium dynamics in such a landscape 
are then of importance: the relaxation time, r r , which is dominated by the curvature of the 
metastable well, and the Kramers escape time from the well, r csc , which involves both the 
characteristics of the well and of the surrounding energy barriers. If the sweep rate, T, of 
the magnetic field is larger than TjT 1 no signature of the metastable well can be seen. If f 
is smaller than t~*, the system relaxes to the equilibrium situation and no magnetization 
plateau is observed. For sweep rates t~* < r < t~ , the system can stay for a while in 
the metastable well and display a M s /3 plateau. Such a picture of the dynamics should be 
supported by a stochastic description of the evolution of a small cluster of spins with a Fokkcr 
Planck equation. In such an approach typically t t = j/loq where ujq is the smallest spin wave 
pulsation in the metastable configuration and 7 is a phenomenological friction coefficient. 
Whatever the origin of the escape from the metastable well, r esc is noticeably larger than 
r r : it involves a factor which depends exponentially on the height, A of the potential energy 
barrier. In the classical regime typically ln(r esc ) is proportional to A/(x"fD), where x is the 
homogeneous susceptibility of the cluster and D the spin diffusion coefficient. In the quantum 
tunnelling regime, which is more probable in the present case, ln(r esc ) oc \J\A\)- More 
theoretical and experimental work is indeed needed to clarify this point. 

In conclusion, we have measured the magnetization process of the S — \ antiferromagnet 



Y. Narumi, K. Katsumata, Z. Honda, J.-C. Domenge, P. Sindzingre, C. Lhuillier, Y. Shimaoka, T. C. Kobayashi and K. Kindo: < 



on the kagome lattice, Cu-titmb, at very low temperatures in both pulsed and steady fields. 
We have found for the first time a M s /3 plateau in a certain range of field sweep rates 
in the pulsed field measurements. A theoretical analysis using exact diagonalization yields, 
Ji = — 19 ± 2K and J2 = 6 ± 2K, for the nearest neighbor and second nearest neighbor 
interactions, respectively. This set of exchange parameters also explains the two-peak feature 
in the magnetic heat capacity observed by Honda et al. [16]. This very strong frustration 
and proximity of the ferromagnetic phase is the origin of a very large density of low lying 
magnetic states. Amongst these states the uud configurations are not very far in energy from 
the absolute ground state and can be selected by an effect of order by disorder. This might be 
an explanation of the transient appearance of a metastable plateau at M s /3 with appropriate 
sweeping rates of the magnetic field. 
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Fig. 1 - The crystal structure of [Cu 3 (titmb)2(OCOCH3)6]-H 2 0. 



Fig. 2 - The magnetization versus applied magnetic field curve of [Cu3(titmb)2(OCOCH3)6]-H20 
measured in increasing and decreasing field runs at T = f 00 mK. Inset: The temperature dependence 
of the magnetic heat capacity of [Cu3(titmb)2(OCOCH3)6]-H20 measured at the designated magnetic 
fields. 



Fig. 3 - The magnetization versus applied magnetic field curve obtained in increasing field runs with 
the designated sweep rates at T = 100 mK. 



Fig. 4 - The magnetization versus applied magnetic field curve obtained in a steady field at the 
designated temperatures. 



Fig. 5 - Dependence of the field range, B2 — Bi, at which the m = | plateau appears on the energy 
per spin, e(m) in this state. Inset: e versus m in the thermodynamic limit. 
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